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Virtual links are generalizations of classical links that can be represented by links embedded 
in a "thickened" surface Ex/, product of a Riemann surface of genus h with an interval. 
In this paper, we show that virtual alternating links and tangles are naturally associated 
with the 1/N 2 expansion of an integral over N x N complex matrices. We suggest that it 
is sufficient to count the equivalence classes of these diagrams modulo ordinary (planar) 
flypes. To test this hypothesis, we use an algorithm coding the corresponding Feynman 
diagrams by means of permutations that generates virtual diagrams up to 6 crossings 
and computes various invariants. Under this hypothesis, we use known results on matrix 
integrals to get the generating functions of virtual alternating tangles of genus 1 to 5 up 
to order 10 (i.e. 10 real crossings). The asymptotic behavior for n large of the numbers 
of links and tangles of genus h and with n crossings is also computed for h = 1, 2, 3 and 
conjectured for general h. 
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0. Introduction 

Virtual knots have been introduced by Kauffman [|l| as an extension of classical knots. They 
may be defined as equivalence classes of 4-valent ( "4-regular" ) diagrams with the ordinary 
under- or over-crossings of knot theory, plus a new type of virtual crossing, depicted with 
a small circle around the intersection, see fig. 1. 




Fig. 1: Ordinary and virtual crossings 




Fig. 2: Ordinary and virtual Reidemeister moves 

Two such diagrams are equivalent if they may be connected by a sequence of gen- 
eralized Reidemeister moves, see fig. 2. An example of a virtual link is provided by the 
following 




(a) (b) ( £ ) 

Fig. 3: A virtual link (a) in the previous notation; (b) as drawn on a genus 
1 surface; (c) alternative representations, see below. 
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Fig. 4: A forbidden Reidemeister move 

In this paper, we also use a different standpoint and notation, closer to graph theory: 
the ordinary crossings are regarded as vertices of a graph, the latter are regarded as rigid, 
i.e. the cyclic order of lines emanating from them cannot be changed, which essentially 
defines a (combinatorial) map, and the virtual crossings are artefacts forced in the planar 
representation by the connections between vertices. In this picture, the meaning of these 
new Reidemeister moves is clear: lines involving virtual crossings may be freely moved 
across the diagram, while keeping their end-points attached to vertices. Also natural in 
this picture is the impossibility for a virtual crossing to pass a line between two ordinary 
crossings, cf fig. 4. 

Virtual knots (or links) may also be thought of as drawn in the vicinity of a connected 
compact orientable Riemann surface E of genus h, (h for "handles"), i.e. embedded into 
the "thickened" surface Ex/, with / an interval. The ordinary over-/under-crossings 
represent the projection of this knot on E, while the virtual ones represent the crossing of 
strands on different faces of E as seen in perspective. See fig. 3(b) for an example. To 
obtain virtual knots, we must consider equivalence classes of such embedded knots modulo 
isotopy in E x /, and modulo orientation-preserving homeomorphisms of E, and addition 
or subtraction of empty handles, see This means that we are interested in virtual 

link diagrams as drawn on "abstract" Riemann surfaces, i.e. independently of the actual 
embedding and without any preferred choice of homology basis (See for example in Fig. 
3(b)-(c) three equivalent representations of the virtual link (a) obtained from one another 
by various modular transformations). As we shall see, this is precisely what Feynman 
diagrams of a matrix integral do for us. 

We shall use for virtual objects the same terminology of knots, links and tangles as 
for classical objects: a link has several connected components, while a knot has only one. 
Tangles, more precisely 4-tangles, have four open ends. 





/ 



\ 



Fig. 5: The flype of a tangle 
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In the same way as alternating knots/links/tangles constitute an important subclass of 
the classical knotted objects, it is natural to define alternating virtual knots/links/tangles: 
they are simply described by diagrams of the previous type, with the condition that along 
any strand, one encounters alternatingly under- and over-crossings, ignoring possible vir- 
tual crossings. Of course, these diagrams have to be divided by the equivalence under 
Reidemeister moves or some combinations thereof. For classical alternating links/knots, 
it is a famous result, conjectured long ago by Tait and finally proved by Menasco and 
Thistlethwaite ||, that it is sufficient to consider reduced diagrams and act on them with 
flypes. Flypes are combinations of Reidemeister moves that preserve the alternating char- 
acter, see fig. 5. 

In a previous work || , we have shown that the integral 

Z N (g,a) = J dMexp-Ntr (aMM ] - ^(MM ] ) 2 ^ . (0.1) 

over N x N complex matrices is well suited for the counting of alternating links and 
tangles: for an appropriate choice of a(g), see below, 2-^- lim/v^oo log Zjy(g, a(g)) is the 
generating function of the numbers of alternating tangle diagrams with n 4-valent crossings, 
and eliminating the equivalence under the flypes just amounts to a "coupling constant 
renormalization" ||, as recalled below. In that way the results of || were reproduced. 

It has been known to physicists since the pioneering work of 't Hooft that the large 
N limit of the previous integral may be organized in a topological way. While the leading 
term corresponds to planar diagrams (in fact, drawn on a sphere), the subdominant terms 
of order N~ 2h of log Z(g) are described by graphs drawn on a Riemann surface of genus 
h. It is thus quite natural to expect that they will be in correspondence with virtual link 
diagrams, (or after differentiation with respect to g, with virtual tangle diagrams) and this 
is what we shall discuss and prove in the following. 

This paper is organized as follows. We first recall (sect. 1) the dictionary between 
link/tangle diagrams and the "Feynman diagrams" of the matrix integral (|0.1|) , and the 



necessary steps to eliminate the redundancy in the counting. First remove spurious di- 
agrams, non-prime or "nugatory" in the knot terminology, by a suitable choice of the 
function a(g). Then we must address the burning question of eliminating the remaining 
redundancies. We argue - but this remains a heuristic argument - that the same (planar) 
flypes as for classical alternating links and tangles are still sufficient to remove the redun- 
dancies of alternating virtual diagrams. We call this the "generalized flype conjecture" and 
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defer until sect. 4 arguments in favor of this conjecture. The flypes are taken into account 
by a redefinition ( "renormalization" ) of the coupling constant: g go(g). In sect. 2 we 
apply these ideas to the case N = 1, which corresponds to the enumeration irrespective of 
the genus; this serves as a sum rule in what follows. In sect. 3 the explicit expressions of 
the first four terms in the large iV expansion of F(g, N) are presented, corresponding to 
genus h = 0, 1, 2 and 3 respectively: while the leading term is well known (and has been 
used in j||) and the second one (h = 1) has been derived by Morris ||, the genus 2 and 3 
contributions computed by Akemann and by Adamietz |9| JT0|1 from the work of [11| may be 



less known. Using the results of previous section and computing the lowest contributions 
of a given genus (Appendix A) gives us enough information to completely determine the 
numbers of all virtual diagrams up to 11 crossings for the links, or 10 for the tangles, and 
up to genus 5. We also give tables of virtual alternating links up to four (real) crossings. 
The rate of growth of the number of virtual diagrams of genus 1, 2 and 3 is also derived 
from these expressions, and a general Ansatz is proposed for generic h, following what is 
regarded as standard lore by physicists. Using the generalized flyping conjecture, we then 
obtain the generating functions of virtual alternating 4-tangles for genus to 3. Their 
asymptotic behaviors for large number of (real) crossing are, up to a larger radius of con- 
vergence, the same as the previous ones, and under the very plausible assumption (proved 
in [JT2JI for classical links) that generic links have no symmetry, we can also estimate the 
asymptotic number of virtual alternating links of a given genus. 

Finally, in sect. 4, we recall that there exists a way to encode the relevant Feynman 
diagrams by means of permutations. This is presumably an old idea; in the present context, 
it seems to be due to Drouffe, it was used in and more recently in related topics by 



combinatorialists ||l~4] - |T5|j . Here we use it to set up an algorithm which is able to build 
all virtual link/tangle diagrams (up to six/five crossings). We have used this to test the 
"generalized flype conjecture" by constructing as many invariants as possible, to make sure 
that the objects that cannot be obtained by flypes from each other are indeed topologically 
distinct. We provide samples of the data thus produced (the full output being accessible 
on the web: [http : //ipnweb . in2p3 . f r/^lptms/membres/pzinn/virtlinks|) , and discuss 
the conclusions one can draw from them. 



1. Matrix integrals and virtual links 

1.1. Feynman rules for matrix integrals 

We first recall the diagrammatic techniques to derive a series expansion in g of the integral 



( p.l|) : see [0|T6[ for a general introduction and || for a discussion in the present context. 
The integration measure in ( p.l| ) is 

dM = Yl dUeMij dSmMij . (1.1) 

l<i,j<N 

We shall be mostly interested in the "free energy" 

F(g,a)= Urn 4log|^ (1.2) 

and its derivatives. The constant a can be absorbed in a rescaling M — > a~^M: 



Z(g,a) = a~ i ^Z(Aj,l) (l-3a) 
F(g,a) = F(-^,l) (1.36) 

but it is useful to keep it. 

Define the "propagator" as the inverse of the quadratic form in ( p.lj ), represented as 
) g > > i = jf^Su5jk and the 4-vertex as the tensor ^vr = gN5 q i5jkSi m o~ n p- This four- 

vertex is to be considered as a rigid crossing, which cannot be flipped and in which the 
cyclic order of the lines cannot be changed. In both the propagator and the 4-vertex, the 
small arrows distinguish the row and column indices of the matrices, while the wide one 
distinguishes M from M'. 

The prescriptions to compute the n-th order of the (/-expansion of F, known as Feyn- 
man rules, are as follows: draw n four-vertices, then draw all the topologically distinct 
connected graphs obtained by joining by propagators the double lines emerging from 
these n 4-vertices, while respecting the orientations, and sum over the matrix indices 
i,j,... = l,...,N. Each graph then comes with a weight g n N# /a 2n , where the power of 
iV will be computed below, and a "symmetry factor" , which is the inverse of the order of 
the group of permutations of the lines and vertices which leave the structure of the graph 
unchanged (see below an alternative characterization of this factor). 

When drawing these Feynman diagrams on a plane, one usually encounters topological 
obstructions which force one to introduce additional crossings (over- or under-, it is imma- 
terial). Alternatively, these diagrams may be drawn on a higher genus Riemann surface 
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as can be seen as follows. The identification of matrix indices by the Kronecker delta's of 
the propagators and 4-vertices leaves us with a number j^F of index loops. By pasting a 
domain homeomorphic to a disk to each such loop, we build a discretized Riemann surface 
with n edges, 2n edges and #F faces. 

Thus Feynman diagrams for F(g, a) may be regarded as discretized orientable Rie- 
mann surfaces E. Their faces are oriented by the small arrows carried by the propagators. 
In addition their edges carry the orientation of the big arrows. According to the argu- 
ment of 't Hooft ]7| and following the rules above, if a diagram has j^V = n vertices, 
hence j^E = 2n edges (propagators), and j^F faces, it carries a power of N equal to 
n — 2n + #F = x#(£) = 2 — 2h, the Euler-Poincare characterics expressed in terms of the 
genus h. 

To summarize, we have obtained a topological expansion (which is an asymptotic 
expansion in 1/N 2 ) 



where is the sum over Feynman diagrams of genus h weigthed as explained above. 
Now that this property has been established, we abandon the double line notation and 
return to more conventional notations for Feynman diagrams: we erase the small arrows 
of matrix indices but retain the big ones that encode the distinction between M and . 

We are now ready to build a dictionary with virtual links: due to the "contraction" 
of M and M* through the propagators, Feynman diagrams of the type just discussed 
are naturally endowed with the properties of alternating virtual link diagrams. Thus 4- 
vertices are in one-to-one correspondence with over /under-crossings \/, while the 

virtual crossings admit the alternative representation: \^ = \/ o ^g^. (Beware! a 
virtual crossing is depicted in the graph theoretic way as an under- or over-crossing in the 
Feynman diagram representation; as mentioned above, it is immaterial to draw it either 
way.) Note that in this representation, it is quite natural that these virtual crossings can 
be freely moved around, thus enforcing the virtual Reidemeister moves. Note also that 
this correspondence gives an operative way to compute the genus on which to draw a 
given virtual link diagram, which may not have been obvious in their original presentation 
(whether or not this is the minimal genus on which one can draw the link itself is a subtle 
matter due to the existence of the "real" Reidemeister moves, and it will be discussed 
again in section 4). 




(1.4) 



h=0 
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From the relations n = #F+2h—2 and j^F > 1, we might have expected contributions 
of genus h to occur first at order n = 2h — 1, for diagrams with a single face. This is indeed 
what happens with diagrams related to hermitian matrix integrals. Here, however, the 
orientation of edges by the big arrows induces an additional constraint: adjacent faces 
have opposite orientations with respect to these arrows. This forbids the possibility to 
have j^F = 1 and as a result, genus h > occurs first at order 

n m m(h) = 2h , (1.5) 

i.e. F^ h \g 7 a) starts at order g 2h , or said otherwise, there is no virtual link of genus h 
with less than 2h (real) crossings. One checks that the bound is saturated by computing 
the coefficient, 

w _ (4fr-l)!! 
* 2h ~ 4/i(2/i + l) • U ' bJ 
and it is also possible to determine without much further effort the next term, i.e. the 
coefficient of g 2h+1 



(,„ (4fc+l)ll ^. (-l)L^J /2h\ 4 *fi-l 



see Appendix A for details. 

1.2. Correlation functions 

We are also interested in the "2p-point functions" G2 P (g,a) := (^tr(MM^) p ), in 



particular 



dg 



and 

G 2 (g, a) = - - ^-F(g, a) = -(l + gG±) , (1.9) 
a oa a 



where use has been made of the homogeneity property ( |1.3| ). This same property implies 
that 

G2 ^' a) = ^(^) (i.io) 

G 2p (g) ■ = G 2p (g, l) 

These functions too admit a graphical representation, with similar Feynman rules and 
Feynman graphs with 2p external lines. These graphs are natural candidates for 2p-tangle 
diagrams. For the four-point function G4, we adopt the following convention of orientation: 
external lines may be extended to a circle surrounding the diagram, and the four lines are 
drawn in the NW, NE, SE and SW directions, with the outcoming arrow (i.e. first crossing 
is over-) on NW and SE external lines. 
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1.3. Nugatory crossings and non-prime diagrams 

The two-point function G2 is in particular useful to dispose of all composite links and 
tangles. Irrelevant "nugatory" crossings and non prime diagrams appear as graphs with a 
subgraph which may be disconnected by cutting transversely two distinct edges. Such a 
subgraph is called a "self-energy" by physicists. To remove all nugatory crossings and non 
prime configurations, i.e. to retain only diagrams with no self-energy, it suffices to choose 
a = a(g) so as to make G^g, cx(g)) = 1 or equivalently in view of ( |1.10| ), / 'ct(g) 2 ) = ct, 
and to plug it into F(g, a) or G2 P (g, a). Then the "connected four-point function with no 
self-energy" defined as 

r(g) = G4g,a(g))-2 (1.11) 

is easily seen to satisfy 

gT(g) = a(g)-l-2g (1.12) 

as a consequence of (|1.9| ). (The appearance of —2 in Ql.llj ) is due to the subtraction of 
disconnected contributions to the 4-tangle by two parallel non-intersecting strands.) 

1.4- Flypes 

We now want to argue that dividing only by the planar flypes suffices to get the equivalence 
classes of virtual alternating tangles and links. This claim is based partially on our intuition 
that other types of moves, such as flypes of higher genus, are not permitted by the structure 
of the thickened surface, and partially on the study of low order virtual links and tangles. 
But our best evidence comes from the analysis, explained in section 4, of several classes of 
invariants applied to links and tangles up to order 6 (six real crossings), which indicates 
that the remaining objects are indeed topologically inequivalent. Still this remains an 
assumption. . . 

For classical (genus 0) tangles, Sundberg and Thistlethwaite have shown how to con- 
struct the generating function of hype-equivalence classes of tangles (g) from the planar 
generating function T^°\g) = lini/v^oo T(g) H- The operations leading from T^°\g) to 
r(°\g) have subsequently been shown by one of us to be simply expressible in terms of 
a "coupling constant renormalization" in the language of physicists, i.e. of a redefinition 
of the expansion variable, determined in a self-consistent way. Let go(g) be the solution of 



Then we recover the result of Sundberg and Thistlethwaite by writing that T^°\g) = 
r(°)((?o (<?))• According to our conjecture, we must generalize this to all genera. If T(g) = 
Sh°=o N~ 2h r( h \g) is the full generating function of virtual prime alternating tangles, then 
we are led to 

r(g) = T(g (g)) . (1.14) 

2. Sum over all genera 

For N = 1 the integral over a single real variable is readily computed. One finds for Z the 
asymptotic expansion 

^-i)i^=f;(i)"^- 

Its logarithm 

F(g, a = 1) | N=1 = g + | g 2 + f g 3 + a$a g 4 + g 5 + g G + m / 

_|_ 14876033 g& _|_ 288018721 g9 _|_ 1227782785 _|_ 142882295557 _)_ Q (^ 12 ) 

(2.2) 

is the sum of for all genera h and will provide a sum rule over the contributions of 
the different genera to be discussed in the next section. By differentiation one gets the 
4-point function 

G 4 (g) = 2 + 10 g + 74 g 2 + 706 g s + 8162 g 4 + 110410 g 5 + 1708394 g 6 (2.3) 
+ 29752066 g 7 + 576037442 g 8 + 12277827850 g 9 + 285764591114 # 10 + 0(g 11 ) 

and after determination of a(g) as explained in sect. 1.3, one gets the connected four-point 
function with no self-energy 

T(g) =2 g + 10 g 2 + 82 g 3 + 898 g 4 + 12018 g 5 + 187626 g 6 (2.4) 
+ 3323682 g 7 + 65607682 g 8 + 1424967394 g 9 + 33736908874 g 10 + <D(g 11 ) . 

Expanding Eq. ( |1.13|) to the required order, we find 

go = g - 2 g 3 - 4 g 4 - 10 g 5 - 30 g 6 - 108 g 7 - 436g 8 - 1890^ 9 - 8588# 10 + 0(g u ) (2.5) 

and therefore 

f(g) =2 g + 10 g 2 + 78 g 3 + 850 g 4 + 11426 g 5 + 179238 g 6 (2.6) 
+ 3187002 g 7 + 63095526 g 8 + 1373767142 g 9 + 3259401885 g 10 + 0(g u ) . 

The three expansions G^g), T(g) and T(g) have the same asymptotic behavior up to 
a multiplicative constant, with their n-th order of the form const. 2 n \Jn + 1 (n + 1)! 
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3. Genus 0, 1, 2 and 3 

We now return to the matrix integral ( p.l|) and its 1/N 2 expansion. Let a 2 (g) be the 
solution of 

a 2 = 1 + 3o(a 2 ) 2 (3.1) 

with a 2 (g) = 1+0 (g). (Its interpretation is that it characterizes the support of the limiting 
distribution of eigenvalues of MMt) Then one finds 

F (0) {g) F (0) {gj 1} = lQga 2 _ _^ (a 2 _ 1)(9 _ fl2) 



n (2n~l)H 
n\(n + 2)! 



n=l 

(As a side-remark, we recall that this is twice the result for the Hermitian matrix integral as 
in (|0.1|) but with tr [% M 2 — f M 4 ) in the exponential). Explicitly, one gets the expansion 

^ (0) (g) = 9 + ! g 2 + 9 g 3 + ^ g 4 + *p g 5 + ^ g 6 + « g 7 

(3.3) 

+ m^m g & _|_ 966654 g 9 + mtmii g 10 + bowsmm g ii + 0(g 12 ) . 
For genus 1, Morris gives [|J 

FW( g ) = -llog (2 - a2)(2 + a2)3 

\HI 24 6 27 

= 1 (3^)- +1 ^ (2n + 2)! _ 

24^ n+1 ^ (n - p)\(n + 2 + p)\ ^ { ' ' 

= ig 2 + m g 3 + mi g 4 + 428 g 5 + ^457 g e + 526I2 g 7 + Mmms ^8 

_|_ 58911256 g 9 _|_ 734641583 _|_ 827733428 fl" + 0(q 12 ) . 

For higher genus, the expressions are more and more complicated. Akemann and 
Adamietz @,|10[ have found that in terms of I\ = 1 — 6ga 2 and Mq = 1 — 2ga 2 



(2) 21aV 69g 2 53g g 



40Jf 640/ 4 2560a 2 / 3 256a 2 / 2 M 

3g 1 3 53 



512a 2 M 3 512a 4 /iM 1024a 4 M 2 15360a 4 / 2 

whence 

F (2) (g) = ?±g 4 + ±ftg 5 + ^g 6 + 46434 ^7 

(3.6) 

+ o 8 + 13798410 g 9 + 10916 r 10282 o 10 + 3328687092 g 11 + 0(g 12 ) 
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and the expression of (g) is too cumbersome to be given here but leads to the expansion 



F (3) (g) = g 6 2504115 g s + 9322668 g 9 + mif*m g 10 + 5345163216 g n + 0(g 12 ) . 

(3.7) 

We check that these F^ h > start at an order in g consistent with ( |1.5| ). Moreover the 
sum of these four first contributions differ from the sum over all genera (|2.2p by terms of 
order g 8 as it should. Using the additional information of (|1.6j ) and (|1.7| ), one may extract 
the first terms of F^ and F^ 

F (4) = 225225 g s + 136865 3 g s + 1495900107 g io + 3023618067 # n + 0(# 12 ) 

(3.8) 

_ 11904165 ^10 _|_ 4304016990 g^- _J_ (3((^ 2 ) 

By differentiating with respect to g, one gets \9i 1) an d G^\g, 1) according to ( |L9| ) 
and ( |1.8| ). One then determines the double expansion in powers of g and 1/N 2 of a(g) so 
as to remove the self-energies, as explained at the end of sect. 1. We don't display the 
corresponding expansion of a(g) as it may be recovered from eqn (|1.12| ) and the expressions 
of r^ 1 ) below. Using the additional data of ( |3.8|) we may provide the g-expansion up to 
order g 10 of T^ h \g) for h = 0, . . . , 5. 

r<°) (g)=g+2 g 2 +6 g 3 +22 g 4 +9l g 5 +408 g 6 + 1938 g 7 +9614 g s +49335 g 9 +260130 g 10 +O(g 11 ) 

(g)=g+8 g 2 +59 g 3 +420 g 4 +2940 g 5 +20384 g 6 _|_i4 479 g 7 +964184 g 8 +6598481 g 9 +45059872 g 10 +O(g 11 ) 
(g) = 17 g 3 +456 g 4 +7728 g 5 + 104762 g 6 _|_i240518 g 7 +13406796 g 8 + 135637190 g 9 + 1305368592 g 10 +O{g^) 

r^l (g) = 1259 g 5 +62072 g 6 _|_i740158 g 7 +36316872 g 8 +627368680 g 9 +9484251920 g 10 +O(g 11 ) (3.9) 

r (4) (g)=200589 g 7 + 14910216 g 8 +600547192 g 9 + 17347802824g 10 +O(c/ 11 ) 

r< 5 ) (g)=54766516 g 9 +5554165536 g 10 +O(g 11 ) 

are the generating functions of the numbers of connected graphs of genus to 3 with no 
self-energy; they are not yet the generating functions of the number of tangle diagrams, 
due to the flype equivalence and possible other redundancies. One may integrate these 
expressions according to ( |1.8|) to obtain the corresponding generating functions F^ 1 ' (g) of 
virtual alternating link diagrams with no self-energy. 

In figures 6-10 we depict the corresponding diagrams of FW, h = 0, 1, 2 up to order 
4: the corresponding diagrams of V are obtained by removing in all possible non equivalent 
ways one vertex, thus opening the link diagram into a tangle. In these figures, we list in 
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Fig. 6: The genus and 1 2-crossing alternating virtual link diagrams in the 
two representations, the Feynman diagrams on the left, the virtual diagrams 
on the right: for each, the inverse of the weight in F is indicated 




Fig. 7 : same for order 3, genus and 1 

parallel the two notations of Feynman diagrams and of links. In the latter, colors have been 
introduced only to distinguish the different connected components. Each link of order n 
(n crossings) comes with an integer, whose inverse gives its weight in F^ h \ Alternatively, 
the number of distinct contributions that this link gives to T^ h ' after removal of one vertex 
equals 2n divided by this integer. 

The first flypes occur at order 3 for genus or genus 1 in T, see fig. 11. The generating 
functions of flype equivalence classes of classical links is 

f(°)(» = T^°\g ) = #+2# 2 +4# 3 +10# 4 +29# 5 +98/+372/+1538/+6755# 9 +30996# 10 +O(# n ) 

(which is the result of 0). Then, according to our assumption, 

f (D (g) = T« (g ) =g + 8g 2 + 57g 3 + 384g 4 + 2512 g 5 

+ 16158 g 6 + 102837 g 7 + 649862 g 8 + 4086137 g 9 + 25597900 g 10 + 0(^ n ) 
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Fig. 8: same for order 4, genus 



f P) (g) = T (2) (g ) = 17g 3 + 456 g 4 + 7626 # 5 + 100910 g 6 + 1155636 g 7 + 11987082 g 

+ 115664638 g 9 + 1056131412 g 10 + 0(g xl ) 
f < 3) (g) = T (3) (g ) = 1259 # 5 + 62072 g 6 + 1727568 g 7 + 35546828 g 8 + 601504150 g 9 



rW (#) = r( 4) (g ) = 200589 g 7 + 14910216 g 8 + 597738946 g 9 + 17103622876 g 10 + 0(g xl ) 
f ( 5 ) (g) = T (5) (g ) = 54766516 g 9 + 5554165536 g 10 + 0(g u ) 



are the generating functions of flype-equivalence classes of virtual tangles with up to 10 
real crossings. For example, the reduction from 59 to 57 of the number of genus 1 tangles 
of order 3 is in accordance with the equivalence of fig. 11. 

The large order behavior of the g-expansions of F^ h > is dominated by the leading 
singularity of a 2 (g), which occurs at g = g c = 1/12, a 2 (g c ) = 2. One finds 



+ 8854470134 g 10 + 0(g 11 ) 



(3.10) 



F^( 9 ,l)^( 9c - 9 f 2 
(g, 1) ~ log (g c — g) . 



(3.11) 



The leading singularity of the expression of comes from the term proportional to i^ -5 , 
as ii vanishes like (g c — g) 1 ^ 2 . Similarly has a pole of order 10 in I\. This is typical 
of what is expected for generic genus WJ\ 



F^(g)^(g c -g) 



5/2(1-/1) 



(3.12) 



Thus one expects 



CO « _ n 5/2(/ t -l)-l 



(3.13) 
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Fig. 10: ibid for genus 2 




Fig. 11: First occurences of flype equivalence in tangles with 3 crossings 
and correspondingly, for T^ h \g) = $^ n 7n <7 n 

7 W « i-„5/2(fc-i) . (3.14) 

The subsequent reductions that we perform to eliminate the redundancies, change the 
value of g c (enlarging it so as to increase the radius of convergence) but do not affect the 
value of the "critical exponents" 5/2, 0, —5/2, . . . Thus, removing the self-energies has the 
effect that the closest singularity is now for g/a 2 (g) = 1/12 which gives [|J g' c = 4/27. 
Similarly, taking care of the flype equivalence increases the radius of convergence of T(g) 
to the value g'£ = (-101 + V21001)/270, see [g, but does not affect the general form ( gTT|) . 
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Finally, note that the function F(g) obtained by integrating V is not exactly the 
generating function of hype-equivalence classes of links (due to the issue of symmetry 
factors) but should have the same asymptotic behavior: intuitively, this reflects the fact 
that the number of link diagrams with a non trivial symmetry factor is subleading and 
does not contribute to the asymptotic behavior of the form ( |3.13| ) (this fact has now been 
proved for classical links [fL2|). 



4. An algorithm to classify virtual alternating links 

In this section we describe the encoding we used to represent virtual alternating link 
diagrams, and the subsequent algorithm that allowed us to generate prime alternating 
virtual links. Due to the factorial growth of their number, we only describe the result up 
at order 4; but we have obtained data up to order 6 in order to check our generalized flype 
conjecture. 

4-1. Alternating link diagrams and permutations 

Our encoding of virtual alternating link diagrams is based on a well-known correspondence 
between bicolored maps and permutations. An alternating link diagram can be equivalently 
described as a (not necessarily planar) map whose vertices have valence 4 and whose faces 
are bicolored, according to the pattern of under /over crossings as one moves around the 
face, see Fig. 12. 



i 




Fig. 12: Bicoloration and labelling of an alternating link diagram. 

Let us now label the edges of the diagram (or of the map): the set of edge labels 
will be called E, and its permutation group S(E). In the implementation, E is chosen 
to be E = {0, 1, . . . , 2n — 1}. It is known that general face- bicolored maps (i.e. duals of 
bipartite maps) are in one-to-one correspondence to pairs of permutations (cr, r) G S(E) 
according to the following recipe: the cycles of a (resp. r) are the labels of the edges in 
their cyclic order as one turns clockwise around white faces (resp. counterclockwise around 
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black faces). □ Define additionally p = a~ x r and p = ar~ l . The cycles of p (or of p) are 
easily seen to be in one-to-one correspondence with vertices of the map. 

Finally a relabelling of the map is a permutation g G S(E) of the labels acting by 
conjugation: 

a' = gag~ 1 r' = grg' 1 (4.1) 

An unlabelled map can therefore be described as a conjugacy class of pairs of permutations. 

Here we require various additional properties of the map, which must be translated 
combinatorially into properties of the permutations: 

1) First and foremost, all vertices must have valence 4. This implies that p (resp. p) 
only has 2-cycles, i.e. is a fixed point-free involution, exchanging edges at overcrossings 
(resp. undercrossings) . Here we decide to focus on p rather than p. The situation at each 
vertex is described on Fig. 13; the figure can be considered as the defining rule to build a 
and r. 

r(a) 



a p(a) 



Fig. 13: Configuration at a vertex. If /3 = p(oe), consistency implies that 
cr(ct) = t(P) and r(ct) = cr(f3) i.e. that p = T~ x a = <j~ 1 t is an involution. 

One can partially fix the freedom on the labels by noting that via conjugations, one 
can reduce p to a given form; for E — {0, 1, . . . , 2n — 1}, we choose 

p(2a)=2a + l p(2a + 1) = 2a a = 0, . . . , n - 1 (4.2) 

Once p is fixed, the data of a alone suffices to describe the alternating link diagram 
since r = ap. Furthermore, all relabellings must commute with p; they form a group 
G = {g G S(E) | gp = pg} which is isomorphic to S n x 2 2 n . 

For example, the labelled diagram of Fig. 12 is such that p is of the form of Eq. ( |4.2j ), 
and we find 



01234567\ /01234567 
'34162705] 7 V43617250 



(4.3) 



1 In terms of the Feynman diagrams of the matrix model, a and r correspond to following 
either of the two lines of an edge in the direction of its (big) arrow. 



17 



or in terms of cycles a = (0 3 6)(1 4 2)(5 7) and r = (0 4 7) (2 5 6)(1 3). 

2) We are interested in connected maps. This amounts to requiring that the action 
on E of the group generated by a and r be transitive. 

3) We mostly focus on diagrams without self-energy. In order to find self-energies (i.e. 
subdiagrams with 2 external legs), we look for pairs of edges (ct,/3) which belong to the 
same cycle of a and to the same cycle of r. Cutting these two edges amounts to composing 
with the transposition (a/3): a' = a o (a/3), r' = r o (a/3). The diagram has no self-energy 
iff for all such pairs (a, 0), the modified diagram corresponding to (a', r') is still connected 
(note that in the planar case the resulting diagram is necessarily disconnected, so that the 
existence of such a pair is enough to discard the diagram). 







i a " 


W(3j 



Fig. 14: Potential self-energies and how to cut them out. 



4) Finally, we want to consider classes of flype-equivalent diagrams. The flype acts 
on a diagram as follows: consider four edges (a,/3,7, 5) in the configuration depicted on 
Fig. 15a), that is o~(a) = /3, 5 and 7 in the same cycle of cr, a and S, (3 and 7 in the same 
cycles of r. Cut the tangle by composing a and r with appropriate cycles, and paste its 
legs together in the way described on Fig. 15b). Proceed only if the resulting subdiagram 
is planar. If it is, then "flip" it by replacing a and r with their inverses inside it, see 
Fig. 15c). Finally, reconnect the tangle to the rest of the diagram, see Fig. 15d). A similar 
operation can be performed by exchanging black and white colors, i.e. a and r in the 
construction above. Together these two types of moves reproduce all possible flypes. 




Fig. 15: Performing a flype via permutations. 
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4-2. Calculation of link invariants 

There are various quantities one may want to compute once a permutation a has been 
produced. They can be of many different types: first, they may be true invariants of 
virtual links, or they may be flype-invariant and therefore conjectured invariants of reduced 
alternating virtual link diagrams, or they may be not invariant at all (but still interesting 
to compute). Secondly, they may be invariants of unoriented or oriented links. In all 
that preceded we have only dealt with unoriented objects; however many useful invariants 
depend on orientation and it is therefore necessary to consider every choice of orientation 
(2 C where c is the number of connected components) of an unoriented object. We now 
list the quantities we have been able to compute, and how to extract them from the 
permutation a: 

(i) The number of crossings n is of course not left invariant by Reidemeister moves, but 
it is preserved by flypes. For reduced alternating diagrams of virtual links it is conjectured 
to be the minimal number of crossings. 

(ii) The genus h of the underlying surface: it is not left invariant by general Rei- 
demeister moves, as Fig. 16 shows (intuitively, after a Reidemeister move a handle may 
become empty so that it must be removed), however it is preserved by flypes, and once 
again, conjectured to be the minimal genus for virtual alternating links. It is given by the 
Euler-Poincare formula: xe = 2 — 2h = j^V — j^E + where j^V = n is the number of 
vertices, #E = In is the number of edges, j^F is the number of faces. If j^a is the number 
of cycles of a i.e. of white faces, and similarly for r, we have #F = j^a + j^r and then we 
conclude that 

fr=l-I(#<7 + #T-n) (4.4) 




Fig. 16: A genus 1 diagram which turns out to be a trivial knot 
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(iii) The number of connected components c is of course an invariant of unoriented 
virtual links. Moving along each connected component on the diagram can be achieved by 
acting with p and p alternatingly; it is easy to check that this implies 

c = \#{P~P) = \#{*- 2 t 2 ) (4.5) 

(iv) The order of symmetry of the diagram: this is not an invariant at all. It is simply 
the order of the group of permutations H that commute with both a and r, that is 

H = {g e G | go = ag} . (4.6) 

This order is a divisor of 2n (this results from the fact that for tangles - see below - this 
group is trivial). 

At low orders one can easily find pairs of flype-equivalent alternating reduced diagrams 
with distinct symmetry factors; this observation is important because it prevents us from 
computing the generating function of the number of prime alternating links in the same 
way as for tangles. 

(v) The set of linking numbers: for an oriented diagram, define the sign e v of a vertex 
v according to Fig. 17, and 

U 3 = ^ l<i,j<c (4.7) 

veVij 

Vij being the set of vertices where components labelled i and j meet. 




Fig. 17: Sign of a vertex. 

The off-diagonal elements £{j, i < j are twice the usual linking numbers between 
components % and j (but they are not necessarily even for virtual diagrams). They are 
invariants of oriented links (up to permutation of the labels of the connected components) . 
Their absolute value does not depend on orientation. 

(vi) The determinant d is an invariant of unoriented links: it is a specialization of the 
usual Alexander polynomial (see below). 
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(vii) The bracket polynomial is defined for an unoriented diagram by a sum over 
"splittings" : 

(L) = J2A a( - s) - b(s) (-A 2 - A' 2 )* 8 ' 1 (4.8) 

s 

where the splitting s is described at each vertex by Fig. 18; a(s) and b(s) are the number 
of vertices of type (a) and (6), and #s is the number of loops thus created.! Alternatively, 
s has a simple description as a permutation: define 

/ .v f a(i) if % overpasses a vertex of type (a) , . _v 

sW = < )■{ - f ■ + r J )>/ (4.9) 

[ t[i) it i overpasses a vertex or type (b) 

for all edges i e E. Note that this induces an orientation of the loops (which is associated 
with the bicoloration of the faces). Then is the number of cycles of s. 

= a ^ + A_ C 

(a) (b) 

Fig. 18: Splitting at a vertex. 

The bracket polynomial is preserved by fiypes; however it is not invariant under Rei- 
demeister move I. It is only up to multiplication by —A that it is an invariant of unoriented 
links. 

One can get rid of this arbitrary power of —A by introducing 

V = (-A)- 3t (L) (4.10) 

where t = Xa<i< c ^ s ^ ne twisting number of the link (which is orientation independent). 
V is (up to a power of A for multi-component links) the Jones polynomial in the variable 
A = x 1 / 4 . 

Furthermore, one can compute colored Jones polyonomials by using cabling, i.e. re- 
placing each string with k parallel strings and then adding extra "twists" to keep constant 
the linking number of each new string with the original one (e.g. keep it zero). We skip the 
details of the implementation; let us simply note that the computation time of the k-th 



2 It might seem surprising that loops that have non-trivial homology are not distinguished 
in Eq. ( |4.8| ); this is because homeomorphisms of E and addition/subtraction of handles do not 
preserve the homology class of loops. 
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cabling roughly grows like c 2 ^ where cik 
in a reasonable amount of time. 



(4fc)! 
(2fc)!(2fc+l)! ' 



so that only k = 2 can be achieved 



(viii) The Alexander polynomials are polynomial invariants of oriented links up to 
a sign and multiplication by a monomial. We refer to [p~8| , pT9| , |T], pO| , [2T[1 for details. The 
(extended, multi-variable) Alexander module is defined by its generators, the edges of the 
diagram, and local linear relations at each vertex, see Fig. 19. They are very simple to 
build in terms of the permutation a once an orientation has been fixed. The th polynomial 
(which vanishes for classical links) is simply the determinant of the matrix of relations. 
Further polynomial invariants are obtained as g.c.d. of minors. 



(1) t h a+(l-st a )b 

(2) bab" 1 




s _1 b t^SL+ (1 - s'H-^b 
b b _1 ab 




Fig. 19: The rules defining (1) the Alexander module and (2) the group of 
a link. 



(ix) Closely related is the group tc of a virtual link which is the generalization to 
virtual links of the fundamental group of the complement of a link; the relations at each 
vertex are described on Fig. 19. It is an invariant (up to isomorphism) of unoriented links 
(the orientation only fixes the presentation). In practice it is not a simple task (and not 
easy to implement by computer) to decide if two groups given by generators and relations 
are isomorphic, and one uses as invariants the number of morphisms of tt into given finite 
groups r. Unfortunately this is only doable for T of small order, which only uncovers a 
small part of the structure of tv. 

4-3. Generalization to tangles 

The appropriate way to consider a (four-legged) tangle is a link with a marked rigid 
vertex. All that has been done in sections 4.1 and 4.2 can therefore be adapted to the case 
of tangles. 

A tangle is now represented by a pair of permutations (a, r) in which the marked 
vertex is encoded just like an ordinary crossing, except the labels must somehow determine 
uniquely which vertex is marked: in the implementation we chose the edges over crossing 
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at the marked vertex to be (2n — 2, 2n — 1). The actual tangle is obtained by removing the 
marked vertex and drawing the external legs in such a way that 2n — 2 is the lower left line 
and 2n — 1 is the upper right line (and therefore, a(2n — 1) is lower right and r{2n — 1) is 
upper left). 

The group of relabellings is now limited to those elements of g leaving the rigid vertex, 
i.e. the corresponding two labels (2n — 2 and 2n — 1), invariant: it is a subgroup G\ of 
G which is isomorphic to S n -i x Z2 n_1 . Other operations on the permutations can be 
performed in the same spirit, that is by keeping the rigid vertex fixed. 

Tangle invariants are generically obtained by pasting an arbitrary given tangle to the 
tangle under consideration and computing link invariants. In practice one obtains only a 
finite number of independent invariants. For example, one gets two Jones polynomials, and 
more generally cik for the /c-cabling, by pasting arbitrary arch configurations to the tangle; 
five th Alexander polynomials obtained by setting equal to two of the four generators 
corresponding to external legs and computing the resulting (n — 1) x (n — 1) determinant 
of the matrix of relations (there are = 6 possibilities but only 5 are independent due 
to a bilinear identity satisfied by the determinants) ; etc. 

4-4- Results 

We have written a program that generates all permutations a G S n up to conjugation by 
elements of G (or G\ for tangles), for n < 6. It then selects connected prime/reduced 
diagrams, creates flype equivalence classes, and finally sorts them according to their in- 
variants and in particular detects undistinguishable non flype-equivalent diagrams. We 
have provided a sample of the output on Fig. 20. In particular, the number of mor- 
phisms of the group 7r into the three groups S3, A4 and A$ is listed, but we have oc- 
casionnaly looked at higher groups. To save space, the variable A of the Jones poly- 
nomial is called a in the tables, and only th extended Alexander polynomials, depend- 
ing on c + 1 variables £q, ■ ■ ■ ,t c _i and s are listed. The rest can be found on the web 
( |http : / / ipnweb . in2p3 . f r/ ^lptms/membres/pzinn/virtlinks|) . 

Up to order 4, we have checked that all flype-equivalence classes are distinguished by 
the invariants described in section 4.2. Therefore the conjecture holds true for links at 
least up to order 4. 

When we go to higher orders, a new problem arises: the difficulty in distinguishing 
links obtained from one another by discrete symmetries. 
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Fig. 20: Table of prime alternating links with 4 crossings. Mirror images are 
indicated only for chiral links. Only th (extended multi-variable) Alexander 
polynomials are listed. 
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Fig. 20 (continued) 

Several discrete operations may be performed within the class of alternating virtual 
link diagrams. First, there is the mirror symmetry, a reflection in a plane orthogonal to the 
plane of the figure. This is a well known operation on classical knots and links, which are 
called achiral or chiral depending on whether the mirror image is equivalent or not to the 
original. The same applies to virtual diagrams. In terms of permutations it corresponds 
to a <-> r. Secondly, there is the simultaneous change of all over- into under-crossings 
and vice versa (which corresponds to a — > r _1 , r — > cr _1 ). For classical links, this is not 
independent of the mirror symmetry, as the composition of the two (a — > cr _1 , r — > r _1 ), 
which is equivalent to a global flip of the diagram, i.e. a rotation around an axis in the 
plane of the figure, yields a link equivalent to the original: this may be seen by gradually 
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overturning the link diagram. For virtual diagrams, this is no longer the case: there is 
an obstruction to this overturning due to the virtual crossings, and the impossibility of 
performing the "forbidden Reidemeister move" of Fig. 4. Accordingly, there are now some 
virtual links that are equivalent to their flip, and some that are not, the latter appearing 
at order 5, see for example Fig. 21. 

In principle, mirror symmetry can be detected by the Jones polynomial since it cor- 
responds to the transformation A — > A -1 . In fact, even the usual Alexander-Conway 
polynomial can distinguish mirror symmetric links in higher genus, since it is no longer 
reciprocal. It is important to notice that unlike classical alternating links, virtual alter- 
nating links do not necessarily saturate the bounds on minimum and maximum degrees 
of their Jones polynomial (maximum degree minus minimum degree in A 4 is less than or 
equal to n — h). For example, there are many virtual alternating links with trivial Jones 
polynomial. Therefore even detection of mirror symmetry can be tricky in higher genus. 
The situation is worse for the flip symmetry since the Jones polynomial (and cabled Jones 
polynomials) cannot distinguish flipped images at all. The group of the link or the higher 
Alexander-Conway polynomials may in some cases distinguish them. 

In practice, already at order 5 there are several diagrams which are not related by 
flypes to their flips but for which we have not found any invariants to distinguish them. 
This is the case for four pairs of links with five crossings, namely those of Fig. 22 and 
23 and their mirror images. At order 6, there are cases of undistinguishable flips and of 
undistinguishable mirror images. Based on the experience at genus 0, we believe that these 
issues are probably difficult to resolve and leave them to future work. 




Fig. 21: A pair of virtual flipped knots, distinguished by their Alexander 
Conway polynomial. 
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Fig. 22: A pair of virtual nipped knots of genus 1, conjectured to be non 
equivalent. 




Fig. 23: A pair of virtual flipped knots of genus 2, conjectured to be non 
equivalent. 

Let us now discuss the case of tangles, which is in fact more important for us since 
it is the objects which we enumerate. Fortunately the rigid vertex destroys any possible 
symmetries and the classification problem becomes easier. Fig. 24 shows a very limited 
sample of our data. 

We have compared the 13010 virtual prime reduced alternating tangle diagrams up 
to order 5. We have performed the flype equivalence and checked that the number of 
tangles of each genus agrees with Eq. (|3.10| ). Our program then allows us to assert that 
all the flype-equivalence classes thus obtained, irrespective of genus and crossing number, 
are distinct (note that in a few cases we had to manually feed the computer groups of 
fairly large order - up to 432 - to make it distinguish the corresponding link groups via 
their morphisms). This we consider a strong argument in favor of our generalized flype 
conjecture. 
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Fig. 24: Table of prime alternating tangles with 2 crossings. * is the type 
of the tangles which encodes how external legs are connected to each other, 
according to: NW connected to t = 1 SE, t = 2 NE, t = 3 SW. 
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Appendix A. Computation of the lowest order term of genus h 

According to the description of the g expansion in terms of permutations of S2 n , the leading 
term with n = 2h reads 

F ^ = (2n- IV 2^1 5 M,{2n}5[ T ], { 2n}5 [(J -i r ] ){ 2™} (A.l) 

where the factor ^r^r comes from the n-th order of the g expansion, and ( 2 n-i)!! takes 
care of the remaining relabeling invariance of the permutations. Also [a] denotes the class 
of the permutation a. One then uses the ortho normalized characters of the S^n symmetric 
group to represent the conditions on a or r by 5[ a ],{a} = J^y.XY {ol)xy {[&]) where the 
sum runs over Young tableaux with 2n boxes, and v a is the number of elements of the 
class a: if a = {l ai 2 Q2 • • •}, v a = (2n)!/ Uji^-j^)- Thus 

(4hVF (h) - v (Ym\ 2 (^L\ 

{ ] 2h ~„k., \m) \{4h)i) 

£ Xy 1 (M)xy 1 ({4M)xr 2 ([r])xy 2 ({4M)xr 3 ([a- 1 r])xn({2 2 n) 



X 



Yi,Y 2 ,Y 3 



The characters xy({2?i}) for the one-cycle class receive contributions only from the hook 

2n — s 



Young tableaux Y s 



for which 



Xv _({2n}) = (-l) 

LfJ 



Y f 

, , , .71—1 



({2"} = (-l)L^J 



f2n-D! 

a. 



s!(2n-s-T 



The summation in ( |A.2|) is thus reduced to a sum over s = 0, • • • ,4h — 1. After some 
algebra, one finds the result (|1.6|) . 



The first subleading term i^ft+i * s gi ven i n a similar way by 



12 n 

= 7^rrTTTT9^T E E 5 Hd2n}5[r],{ P) 2n-p}5[ CT -l T ], { 2™} (A.3) 



2n-l)!! 2 n n\ 
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with now n = 2h + 1, since F = 3 implies that either j^a = 1, j^r = 2 or the converse. 
The same method as above applies, again only hook Young tableaux contribute, and the 
only additional piece of information required is 

f (-l) s ifs + l<p<2n-s-l, 
X Ya ({2n-p,p}) = I (-l) s+1 it2n-s<p<s 
1 otherwise 

as we learn from the Murnaghan-Natayama formula. A little algebra then leads to the 
result (fL7D . 
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